Planck2013 results support the simplest cyclic models 
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We show that results from the Planck satelhte reported in 2013 are consistent with the simplest 
cyclic models for natural parameter ranges {i.e., order unity dimensionless coefficients), assuming 
the standard entropic mechanism for generating curvature perturbations. With improved precision, 
forthcoming results from Planck and other experiments should be able to test the remaining pa- 
rameter range and confirm or refute the core predictions - i.e., no observable primordial B-mode 
polarization and detectable local non-gaussianity. A new prediction, given the Planck2013 con- 
straints on the bispectrum, is a sharp constraint on the local trispectrum parameter Qnl', namely, 
the simplest models predict it is negative, with qnl ^ —1700. 



The recent results by the Planck satellite have been 
presented as supporting the predictions of some of the 
earliest and simplest inflationary models of the early uni- 
verse [1-3], which can all be reformulated as single-field 
models with a plateau type potential [4]. However, these 
models have a severe initial conditions problem and typ- 
ically lead to eternal inflation, so that, in claiming to 
match the 'predictions' of these models, one is implicitly 
making a number of unstated and poorly-understood as- 
sumptions about measures and the multiverse [4]. This 
is not an insignificant point: with even moderately dif- 
ferent initial conditions and/or with some of the simplest 
measures (such as physical volume weighting) the predic- 
tions of these simplest inflationary models are in fact in 
complete disagreement with the Planck2013 results. Fi- 
nally, within the context of the inflationary paradigm, 
the observationally favored potentials are highly unlikely 
since they are less generic, require more tuning, inflate 
for a smaller range of field values, and produce expo- 
nentially less inflation than the power-law potentials [4] . 
Thus, looking for alternatives to inflation remains as well 
motivated as ever. 

In particular, the cyclic theory of the universe [5] does 
not appear to suffer from these same problems. In the 
present paper, we demonstrate that the simplest cyclic 
models can match all of the cosmological parameters fa- 
vored by Planck. These models are characterized by sim- 
ple potentials; by natural parameter ranges {i.e., with di- 
mensionless quantities of order unity) ; and they assume 
the standard entropic mechanism for generating curva- 
ture perturbations [6]. At the same time, forthcoming 
results from Planck and other experiments will reach the 
point where they should be able to explore the natu- 
ral ranges and confirm or refute the core predictions - 
namely, no observable primordial B-mode polarization 
and detectable local non-gaussianity in the bispectrum. 
Furthermore, given the Planck2013 constraints on the 
bispectrum, we present a new test, namely a sharper pre- 
diction for the local trispectrum parameter gNL- 



Our conclusions differ markedly from those reported 
by the Planck2013 collaboration who describe the ekpy- 
rotic/cyclic scenarios as being under "severe pressure" 
[3] as a result of the new data. We note that the basis 
of their summary is an analysis of a particular form of 
scalar field potential (exponential); a particular mecha- 
nism of generating curvature perturbations (entropic, as 
in this paper); and, within that mechanism, a particular 
choice of conversion from entropic to curvature perturba- 
tions occurring during the ekpyrotic smoothing phase [7] , 
which, in the cyclic model, is far less likely than the alter- 
native, a conversion after the smoothing phase is over - 
see [8]. A fair analogue would be to judge the inflationary 
scenario on the basis of how well an exponential inflaton 
potential fits the Planck2013 data: the answer would be 
it is ruled out at more than the 99.7% confidence level [2]. 
Obviously, it would be wrong to rule out a whole scenario 
on the basis of a single type of potential. Furthermore, 
in their detailed discussion, the Planck2013 collaboration 
reports that the same data fits the cyclic predictions for 
a range of parameters just by shifting the conversion to 
the kinetic phase occurring after the smoothing phase 
and still keeping to the strict exponential form [3]. This 
Planck2013 result is consistent with our analysis below, 
which assesses the situation introducing a more model- 
independent parameterization of the ekpyrotic potential 
that enables us to quantify the fit in an unbiased way. At 
the same time, we note that a new prediction regarding 
the trispectrum resulting from the Planck measurements 
in the context of the entropic mechanism. More gen- 
erally, given the serious and exciting scientific issues at 
stake, such as the origin and evolution of the universe, 
the nature of the cosmic singularities, and the existence 
of a multiverse versus universe, we hope data analysts 
will take some care not to reach conclusions about a cos- 
mological scenario prematurely. 

Cyclic Model. The cyclic model [9, 10] is based on the 
idea that the big bang is a big bounce from a phase of 
contraction to a phase of expansion; that the key events 



that set the large-scale smoothness, flatness and density 
perturbations occurred before the bounce during a phase 
of dark energy domination (similar to today) and ultra- 
slow (ekpyrotic) contraction with equation of state w ^ 
1 [11] (with no high energy inflation after the bounce); 
and that the bounces regularly repeat such that phases 
of expansion and contraction alternate. Unlike earlier 
oscillatory models, there is a strong asymmetry between 
phases, with each cycle undergoing vastly more expansion 
than contraction. Thus, during each cycle the scale factor 
a grows by a huge factor. Nevertheless, local quantities, 
such as the Hubble parameter H = a/a, temperature, 
and density return to the same value after each cycle, 
and it is in this sense that the models are cyclic. 

The ekpyrotic phase can be generated by a scalar field 
a evolving down a steep and negative potential 



during the ekpyrotic phase as 
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where Vq is a constant, e = |(1 + w) is in general a 
slowly- varying function of time and the ellipsis stands for 
additional terms that parameterize the potential in terms 
of a second scalar field s that we will discuss below. 

The big bounce has been modeled in two inequivalent 
ways so far: in the first class of models the bounce is 
classically non-singular in the sense that the scale factor 
a reverses from contraction to expansion before reaching 
zero. Such models typically rely on higher-derivative ki- 
netic terms [12-17]. In the second class of models, which 
were originally motivated by the braneworld picture sug- 
gested by string theory, the bounce corresponds to the 
collision of orbifold branes. In this case, the evolution is 
classically singular (a reaches zero), but a semi-classical 
treatment indicates that this singularity is resolved at 
the quantum level [18]. In either case, the bounce is an 
analytic and unitary event so that all physical quantities 
vary continuously across the bounce, which is all that 
needs to be known to make predictions. 

Notably, the Hubble rate is largest right at the begin- 
ning of the radiation phase after the bounce, and is signif- 
icantly smaller during the 'smoothing' dark energy and 
ekpyrotic phases. Thus, quantum fluctuations during the 
smoothing phases are always small, and consequently no 
run-away behavior, such as the one leading to eternal 
inflation, occurs in these models [19, 20]. This feature 
makes it possible to avoid the unpredictability problem 
of inflation. 

Density perturbations. The best-understood way in 
which nearly scale-invariant curvature perturbations can 
be generated during a contracting phase of the universe 
is via the entropic mechanism [6, 21]. In this scenario, 
there is a second scalar field s, which develops nearly 
scale-invariant entropy (or isocurvature) perturbations, 
and these entropy perturbations are subsequently con- 
verted into curvature perturbations during the kinetic 
phase just before the bounce. 

We will parameterize the potential for both scalar fields 
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where K2,K3,K4 are constants. A potential of this form 
can for example be constructed by having two scalar 
fields 01 , 02 with separate ekpyrotic potentials V = 
_yj^gci0i „ Y2gC2ci>2 g^j^(^ then performing a rotation in 
field space into the ekpyrotic direction a (defined to 
point tangentially to the background trajectory, with 
& — (01 + 92?'^'^) and the transverse direction s. In 
this case, the K2.3 4 parameters are 



'«2 = 1, K3 = 2v2 



{A-cD 



\c\c2\ 



K4 = 4- 



(cf 



c\c\{c\ + ^^' 



(3) 

with 1/e = 2/cf -I- 2/c|. These expressions allow us to 
estimate natural parameter ranges: we may for example 
expect K2 to be close to 1, say within 10%. One may 
also expect that in a fundamental physics context, the 
values of c\ and ci that one can obtain are not wildly 
different. In this case the potential is roughly symmetric 
across s = 0, and as we will see below, this case is the one 
of most interest. Then we may expect K3 G [—1,1] and 
K4 K, 4. The ekpyrotic phase ends when the steeply falling 
potential bottoms out at Vek-end, a short time before the 
bounce. We are now ready to turn to a calculation of the 
linear fluctuations in these models. 

There are two gauge-invariant scalar fluctuations to 
consider: the curvature fluctuation (, deflned as a per- 
turbation to the scale factor in a flat FLRW universe, 



ds^ 



-dt^ 



a^{t)e^^^'-''^Ux'dx, 



(4) 



and the entropy fluctuation 5s. During the ekpyrotic 
phase, the curvature perturbations obtain a very blue 
spectrum with a small amplitude, and hence can be ne- 
glected [22]. However, the entropy perturbations are of 
more interest. They obey the equation of motion 



5s + ■iH5s + V.ss = 0. 



(5) 



A standard calculation then shows that this equation can 
be solved in terms of Hankel functions, and the bound- 
ary condition that in the far past the solution should 
approximate the Minkowski vacuum flxes the solution to 
be given by 5s = ^~^ HI, {—kt), where the Hankel func- 
tion has index 1^ = {1 + 2k,2-^ - ^ + §^2^)^/2, Here 
Af denotes the number of e-folds left before the end of 
the ekpyrotic phase {dJ^ = d{lnaH)) and terms of order 
1/e^ have been neglected. At late times {—kt) — > we 
obtain 5s « -j^. — r— , implying that at the end of the 

ekpyrotic phase, the entropy perturbation is given by 
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Hence, the spectral index is given by 
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where the last approximation is valid as long as K2 is close 
to 1, i.e. as long as |k2 — 1| ^ 1- Thus, the spectrum 
typically deviates from exact scale-invariance by a few 
percent, either to the blue or to the red. The term in k,2 
can by itself render the spectrum red or blue. The 2/e 
contribution shifts the spectrum to the blue, while the 
e.A^/e term tends to shift it to the red, as towards the end 
of the ekpyrotic phase e can be expected to decrease while 
Af decreases by definition. Thus, the entropic mechanism 
does not yield a clear preference for either a red or a blue 
tilt. However, for natural values of the parameters, the 
spectrum can be in perfect agreement with the measured 
Planck2013 value n^ ^ 0.9603 ± 0.0073 [1], for example 
by having K2 = 1.06 and constant e = 50. (Recall that 
e is typically of order 1/M in inflationary models and of 
order M in ekpyrotic/cyclic models where A/" ^ 50 for the 
modes of observable interest in both cases.) 

Conversion of the entropy perturbations into curva- 
ture occurs gcnerically when the trajectory in scalar field 
space bends (for a concrete realization in supergravity, 
see [23, 24]). The result of this bending is that a (lin- 
ear) curvature perturbation C,l develops, with a spectrum 
identical to that of the entropy perturbation 5s and with 
an amplitude C,l ~ \5s{tek-end)- Thus, the amplitude of 
the resulting curvature perturbations is given by 
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so that the measured variance can be achieved for ekpy- 
rotic potentials bottoming out at Vek-end around the 
grand unified scale [25] . 

Non-Gaussianities. Aside from improving measure- 
ments of the power spectrum, a great advance in 
Planck2013 is a substantially better bound on non- 
gaussianity [3] . Since any non-gaussian signals are seen as 
clear indicators for something other than simple infiation- 
ary models, the results rule out a large range of complex 
inflationary models and other alternatives. However, we 
show below that the current limits fit simple cyclic mod- 
els for a substantial natural range of parameters. 

For cyclic models, one intuitively expects significant 
departures from gaussianity compared to inflation be- 
cause the requisite potential is steeper and thus self- 
interactions of the scalar fields are guaranteed [7, 26- 
28]. But, as we will review below, this does not imply 
that every non-gaussian estimator must be large. In cur- 
rent models, the kinetic terms for the scalar fields are 
of canonical form. Hence, the non-gaussian corrections 
that are generated are of the local form, and they can be 
analyzed by studying the classical equations of motion 
to higher orders in perturbation theory. It is useful to 
define the following expansion for the curvature pertur- 
bation (on uniform energy density hypersurfaces) 
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where the numerical coefficients are conventional. The 
parameters Jml and qml then describe the deviations 



from gaussianity of the bispectrum and trispectrum re- 
spectively. 

In order to calculate these parameters, we first have 
to work out the entropy perturbations to third order in 
perturbation theory. A straightforward, but lengthy cal- 
culation yields [29] 

,,..„+ -^,4 + .(-4 -^,,4, (10, 

where terms of order 1/e have been neglected. In order 
to calculate how this entropy perturbation gets converted 
into a non-linear curvature perturbation, it is convenient 
to use the formula [301 
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which is valid to all orders and where an overbar denotes 
a background quantity. One can use this simple formula 
to derive an analytical estimate of the non-linear curva- 
ture perturbation. However, a more precise procedure 
is to integrate the above formula numerically. This was 
done in [29, 31, 32] for several different implementations 
of the bending of the background trajectory. There it 
was found that the results are most robust when the con- 
version process is relatively smooth (i.e. lasting on the 
order of one Hubble time - shorter conversion times tend 
to amplify the magnitudes of the non-linearities some- 
what), and we will focus on this case. The resulting 
non-gaussianity parameters can be well fitted by the ap- 
proximate formulae 
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which also agree well with the analytical estimates. 

Along with the spectral tilt in Eq. (7), these two ex- 
pressions comprise the key predictions of the scalar den- 
sity perturbations (generated via the entropic mecha- 
nism) to compare with Planck2013 and future observa- 
tions. In addition, ekpyrotic/cyclic models predict no 
observable tensor modes on large scales [33, 34]. 

The expressions here differ in form from previous 
work because we have introduced the three dimension- 
less model-independent parameters k^ to judge the natu- 
ral range expected in this cyclic model. The essential re- 
sult remains: the model predicts substantially more non- 
gaussianity than simple infiationary models. For exam- 
ple, whereas the total f^L for simple inflation models is 
less than one [35] , we see that the formula for /at^ in the 
cyclic model contains a term ±5, which arises from the 
conversion process. Its sign depends on the details of the 
conversion process {e.g. the trajectory bending to the left 
rather than to the right in fleld space). The second term 
in Eq. (12) is also substantially greater than one. How- 
ever, allowing for either possibility {i.e. plus or minus 



sign), we find that the Planck bounds of f^L ~ 2.7 ± 5.8 
at the 1(7 level can be matched respectively by 



-0.76 < K3 < 0.33 or 0.18 < K3 < 1-27 
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where we have used again a typical value of e = 50. These 
allowed ranges for K3 at the la level are still of 0(1), 
and thus, contrary to the claim in [3], do not constitute 
"severe pressure" on cyclic models. The fit remains good 
over a substantial range of e. 

In fact, fine-tuning is not required to achieve any value 
in the entire range of /nl allowed by Planck2013, includ- 
ing fisiL near zero (as expected for simple inflation mod- 
els). Hence, observing I/tvlI measurably greater than 
one would fit the simplest cyclic picture and be incon- 
sistent with simple infiation models; but observing I/nlI 
less than one would, taken by itself, not be conclusive. 

Fortunately, there is a another test to consider: the 
trispectrum. In addition to the parameter g^^ that we 
are about to discuss, there exists a second trispectrum 
parameter, denoted r^vL, which for models where the 
perturbations originate from a single perturbation mode 
(such as the entropy perturbation here) takes the form 
(see e.g. [36]): 
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Thus, given the bounds on Jnl that we just described, a 
consistency check for ekpyrotic models is that tml must 
be smaller than about 100 at the \(j level. The con- 
straints from the Planck satellite for this parameter are 
Tnl < 2800 (95% confidence level), and thus these are 
consistent. 

Less well-constrained, but probably also more inter- 
esting, is the intrinsic local trispectrum parameter gNL- 
The strong bounds on /nl imply that K3 cannot be much 
greater than unity, and in turn this implies that K4 « 4 
(see the discussion below Eq. (3)). Using these relations, 
we can simplify the fitting formula for the trispectrum, 
so that it ends up depending solely on the equation of 
state parameter e. 



9NL 



-35 e. 



(16) 



This remarkably simple formula constitutes a key test 
for the simple cyclic models described here. The percent 
level deviation of the spectral tilt from exact scale invari- 
ance is natural if e is equal to or larger than about 50, cf. 
Eq. (7) - this consideration immediately leads to an ap- 
proximate upper bound on gj\jL of about g^^ ^ —1700. 
A larger equation of state parameter e then leads to a 
correspondingly larger magnitude of g^L- 

What is interesting is that the entropic mechanism 
makes a definite prediction for the sign of gNL , namely 
that it should be negative. We note that a larger e allows 
a smaller range of K3 values to be consistent with the 
current bounds on fNL {cf Eq. (12)), so that one may 
conclude that values of gNL closer to the upper bound. 



— 1700, are more natural than substantially larger neg- 
ative values. The robustness of the result is due to the 
large contribution from the K3 , K4-independent term in 
Eq. (13). Although current bounds on gNL are still weak 
(they are of order ±10^ [37, 38]), this prediction is inter- 
esting for future observations, as the predicted values are 
much larger in magnitude than those of simple inflation- 
ary models. 

Discussion: The data recently released by the Planck 
satellite team is extraordinarily precise. It shows that 
the primordial perturbations were nearly scale- invariant, 
adiabatic and gaussian. These flndings have immediately 
been hailed as confirming the predictions of some of the 
simplest models of inflation. However, as emphasized in 
Ref. [4], at the same time the data favors certain models, 
it disfavors the paradigm. The particular models favored 
by the data are a set that is exponentially unlikely to 
have produced our universe compared to more generic 
power-law potentials and only produce inflation in cases 
where the universe is already surprisingly smooth over 
many Hubble volumes before inflation begins. 

By contrast, we have shown in this paper that a rad- 
ically different class of models without these problems, 
namely cyclic models containing an ekpyrotic phase and 
using the standard entropic mechanism [6], match the 
Planck findings very well: the Planck central values of 
cosmic microwave parameters, as well as the entire la 
ranges , can be reproduced with natural 0(1) ranges of 
the parameters. At the same time, by constraining the 
spectral tilt and fNL to a significant degree, the cur- 
rent results imply that a new (falsifiable!) prediction for 
cyclic models is that the intrinsic trispectrum parameter 
gNL should be negative in sign and of order 10"^ — 10^ in 
magnitude. (As a cautionary note, we point out that the 
method of generating curvature perturbations in cyclic 
models is, to a large degree, separate from the scenario 
overall, as in the case of inflation. Other approaches than 
the one considered here would produce different predic- 
tions for non-gaussianity.) Whereas the non-gaussianity 
is somewhat model-dependent, a model-independent, fal- 
sifiable prediction of cyclic models is that there are no de- 
tectable primordial gravitational waves on the scales of 
cosmological interest. Thus, near-future measurements 
will provide clear-cut tests of these simple cyclic models 
of the universe. 
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